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Abstract 

A  procedure  to  predict  the  state  of  charge  of  a  lithium  ion  cell  using  experimental  data  (cell  potential  versus  time),  as  it  becomes  available,  is 
presented.  The  procedure  is  based  on  the  physics  of  the  system  and  provides  a  realistic  estimate  of  the  state  of  charge  of  the  cell  as  a  function  of 
time,  for  a  given  set  of  properties  of  the  electrodes.  An  electrochemical  cell  model  is  used  to  obtain  an  extended  Kalman  filter  (EKF)  for  estimating 
the  state  of  charge  (SOC)  of  a  lithium  ion  cell  in  which  the  negative  electrode  is  the  limiting  electrode.  The  method  could  also  be  used  for  a  cell  in 
which  the  positive  electrode  is  limiting. 
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1.  Introduction 

The  phrase  state  of  charge  (SOC)  is  routinely  used  to  describe 
the  performance  of  a  lithium  ion  cell  [1].  Numerous  approaches 
have  been  presented  in  refs.  [2-28]  to  monitor  the  SOC  of  a 
cell;  most  of  these  are  for  lead-acid  batteries.  These  include, 
but  are  not  limited  to,  following  a  physical  property  like  the 
internal  resistance  of  the  cell  or  the  electrolyte  density  [3,4], 
measurement  of  the  cell  impedance  [5-8],  optical  methods  [9], 
eddy  current  methods  [10],  equivalent  circuit  analysis  of  the 
charge/discharge  and/or  impedance  curves  [11-18,22],  tech¬ 
niques  employing  a  Kalman  filter  [19-21],  techniques  based  on 
fuzzy  logic  [23-25]  and  coulomb-counting  [26,27].  The  model 
by  Tenno  et  al.  for  the  VRLA  battery  [2]  is  valid  only  under 
overcharge.  Piller  et  al.  [28]  provide  a  good  overview  of  most 
of  these  methods.  Huet  [5]  provides  an  extensive  review  of 
impedance  measurements  for  determination  of  the  SOC  for  lead- 
acid  and  nickel-cadmium  cells.  Verbrugge  et  al.  [14-17]  discuss 
the  utility  of  circuit  analog  models  in  the  recursive  estimation 
of  parameters  to  predict  battery  performance  in  the  context  of 
electric  vehicles.  Fuzzy-logic  based  models  that  can  be  readily 
implemented  in  hardware  are  discussed  in  the  work  by  Salkind 
et  al.  [23].  Hansen  and  Wang  [12]  present  a  pattern  recognition 
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algorithm  that  reduces  the  computational  effort  for  circuit-based 
models  developed  using  ADVISOR  [13]  and  makes  them  suit¬ 
able  for  real-time  embedded  system  applications.  Nanjudaiah 
and  Koch  [26]  developed  a  cycle-life  sensor  for  a  Li/Li+  half 
cell  by  inserting  a  Pt  microelectrode  inside  the  cell.  Martinet  et 
al.  [18]  present  a  correlation  between  the  electrochemical  noise 
in  a  battery  system  and  its  state  of  charge. 

In  methods  employing  circuit  analogy,  the  parameter  that  is 
tracked  to  estimate  the  state  of  charge  is  most  often  empirical 
with  little  physical  significance.  The  coulomb-counting  proce¬ 
dure  is  not  reliable  since  not  all  the  current  supplied  goes  to 
charging  the  cell  and  the  charge  lost  to  any  undesired  side 
reaction  [30]  is  not  included.  Similarly,  the  fuzzy  logic  based 
methods  rely  on  the  training  data  supplied.  Hence  the  results 
from  these  models  are  of  limited  value  in  predicting  the  SOC  or 
SOH  of  a  battery.  Consequently,  there  arises  a  need  for  state  of 
charge  estimation  based  on  the  physics  of  the  cell.  Unfortunately, 
the  computational  time  required  to  solve  rigorous  battery  models 
[29]  repeatedly  practically  eliminates  the  utility  of  physics-based 
models  that  are  used  to  characterize  the  SOC  of  the  cell.  Inter¬ 
mediate  models  (i.e.,  extension  of  circuit  analog  models)  were 
presented  by  Plett  [19-21],  but  these  are  too  empirical  to  pro¬ 
vide  any  insight  into  the  system.  Some  of  our  recent  publications 
[31-33]  present  a  simple,  computationally  efficient,  ordinary 
differential  equation  (ODE)  model  of  the  cell  using  the  FOR¬ 
TRAN  subroutine  DDASL  [34],  that  is  capable  of  predicting 
well,  the  charge/discharge  curves  of  a  lithium  ion  cell,  up  to 
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Nomenclature 

cij  specific  surface  area  of  the  electrode  (m2  m-3) 

c\j  concentration  of  lithium  in  the  solid  phase 

(molm-3) 

co  concentration  of  the  electrolyte  (molm~3) 
c?vf  average  concentration  of  lithium  in  the  solid  phase 

(mol  m  3) 

cmax  maximum  concentration  of  lithium  in  the  solid 

’  J 

phase  (molm-3) 

Ci  ■  concentration  of  lithium  at  the  surface  of  the 

’  J 

sphere  (mol  m-3) 

c?  :  initial  concentration  of  lithium  in  the  solid  phase 

(mol  m  3) 

D\  j  diffusion  coefficient  of  lithium  in  the  solid  phase 
inside  electrode  (m2  s_1) 
e  weighted  sum- squared  error 

/  forcing  function  matrix 

F  Faraday’s  constant  (96,487  C/equivalent) 

F  linearized  forcing  function  matrix 

G(t )  matrix  of  process  noise  coefficients 

h  function  of  the  estimated  states 

H  linearized  system  matrix 

Jj  flux  entering  electrode  j  (A) 

kj  rate  constant  for  the  electrochemical  reaction  at 
the  surface  (m2,5  mol-0,5  s-1) 

K  filter  gain 

n  number  of  electrons  transferred  (=1) 

P  covariance  matrix 

Q  Variance  of  the  process  noise 

r  radial  coordinate  (m) 

Rj  radius  of  the  particle  in  phase  (j  =  F  p)  (m) 
74+i  variance  of  the  measurement  noise 

t  time  (s) 

u(t)  input  function 

Uj  open  circuit  potential  at  the  electrode  . 

x  vector  of  the  state  variables 

Vceii  cell  potential  (V) 

W  weight  matrix 

y  vector  of  the  predicted  values  of  the  output  vari¬ 

ables 

y  vector  of  the  experimental  values  of  the  output 

variables 

01  j  solid  phase  potential  at  the  electrode  (V) 

Oj  utility  coefficient  of  electrode 

6n  state  of  charge  based  on  the  negative  electrode 


and  including  the  1  C  rate.  In  this  work,  our  model  is  used  to 
design  an  extended  Kalman  filter  (EKF)  to  predict  the  state 
of  charge  of  a  lithium  ion  cell  as  it  is  charged  or  discharged. 
The  model  is  regressed  with  synthetic  data  (cell  potential  as  a 
function  of  time)  to  provide  estimates  of  the  solid  phase  con¬ 
centrations  as  functions  of  time  and  hence,  the  state  of  charge  of 
the  cell  can  be  monitored  online.  State  of  health  (SOH)  on  the 


other  hand  provides  information  on  aging  of  the  cell.  Cycling 
performance  of  the  cell  is  not  considered  in  this  paper  and  hence 
SOH  calculations  are  of  little  relevance  in  this  context  [20]. 
However,  it  must  be  mentioned  that  an  extension  of  the  present 
technique  (currently  being  pursued  in  our  laboratory)  to  study 
cycling  behavior  will  yield  the  SOC  as  well  as  the  SOH  of  the 
cell. 


2.  The  Kalman  filter 


The  term  “filter”  refers  to  a  weighted  least  squares  algorithm 
that  can  be  used  to  minimize  the  effects  of  measurement  noise 
[45].  A  brief  review  of  linear  filters  is  provided  in  this  section, 
followed  by  a  direct  analogy  for  the  nonlinear  case,  which  is  of 
interest  in  the  context  of  an  electrochemical  cell. 

Let  us  say  we  have  an  approximate  model  of  a  system: 

y\(t)  =  a\x\(t)  +  b\x2(t) 
y2(t)  —  a2xi(t)  +  b2x2(t) 

where  y\  and  y 2  are  equivalents  (in  the  model)  of  the  variables 
that  can  be  directly  measured  in  an  experiment,  x\  and  X2  are 
the  state  variables  (in  the  model)  to  be  estimated  as  a  function  of 
time,  for  given  values  of  the  constants  at  and  b[.  The  objective 
is  to  estimate  and  update  the  values  of  x\  and  X2  as  we  progress 
in  time. 

The  equation  above  can  be  written  in  the  state- space  form 
[45]  as  follows: 


y  =  H  x 

where  the  output  vector  is: 


yi 

y2 


and  the  vector  of  states  is: 


H  is  the  system  matrix  (or  the  matrix  of  coefficients): 


b\ 

b2 


Let  y*  and  y |  be  the  experimental  values  corresponding  to  the 
model  predictions  y  1  and  y 2,  respectively.  From  an  experiment, 
we  are  able  to  measure  and  y|  at  various  time  steps,  for 
example,  4  and  4+1.  Let  us  denote  these  values  as  y\k  and 
y\  k  at  t  =  4,  and  y\  k+l  and  y|  £+1  at  t  -  4+1  •  The  corresponding 
values  in  the  model  are  y\^  and  y2f  for  t  -  4  and  y \ ^+1  and  y2,k+\ 
for  t  =  4+i. 

In  a  set  of  experimental  data,  all  the  data  points  need  not  be 
equally  weighted — some  allowance  must  be  provided  to  account 
for  the  deviations  due  to  experimental  error.  Let  the  weight  asso¬ 
ciated  with  a  data  point  at  time  t  =  4  be  Wkk  (he.,  the  weight  for 
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y*  k2Ltt  =  4  will  be  Then  the  weighted  sum-squared  resid¬ 

ual  is  given  by: 

ek  =  y(yk-ypTwk(yk-yp]  (6) 

where 


is  the  filter  gain  and 


Pk+l  =  ( HjWkHk  +  Hj+lWk+lHk+l) 
is  the  covariance  matrix. 


(17) 


assuming  that  there  is  no  correlation  between  the  errors.  Note 
that  if  7’  is  the  number  of  state  variables,  then  the  residual  is 
given  by: 


3.  Extension  to  the  nonlinear  case 

For  the  nonlinear  case,  the  model  equations  are  typically  of 
the  form: 

x(t)  =  f(x(t ),  ii(0,  0  +  G(t)W(t)  (18) 


n  =  - y*P2  (8) 

j 

Our  objective  then,  is  to  estimate  the  state  variables  (x\tk  and 
X2,k)  such  that  the  error  between  the  predictions  made  (y 14  and 
y2,k)  using  the  estimates  and  the  experimental  measurements  is 
minimal:  i.e.,  minimize  the  expression  for  the  sum  squared  error 
given  by  Eq.  (8),  with  respect  to  the  state  variables.  To  achieve 
this,  we  find  the  first  derivatives  of  the  function  to  be  minimized 
(ek)  with  respect  to  each  of  the  state  variables  (x\ 4  and  X24)  and 
set  them  equal  to  zero: 


1  v-^  dyjk 

v»«  = 


WjAyj,k-ylk)  = 0  for*  =  1,2  (9) 


7=1 


Rewriting  Eq.  (9)  in  the  state-space  form,  we  have: 

Vxkek  =  HjWk(yk-y*k)  =  0  (10) 


Since  J4  =  //£  x  Xk  Eq.  (10)  can  be  rewritten  as  follows: 

Vxkek  =  Hj Wk(Hk  ■xk-y*k)  =  0  (11) 


Eq.  (11)  yields: 

xk  =  ( HjWkHkf1  Hjwky*k  (12) 

We  now  have  Xk  in  terms  of  y|.  The  next  step  is  to  find  Xk+i  in 
terms  of  Xk  and  y^+1 .  At  time  4+1  we  have  data  points  from  4 
and  4+i .  Hence  using  similar  arguments  [45]  used  to  derive  (12) 
we  have: 

•*4+1  =  (■ HjwkHk  +  Hj+lWk+ 1  Hk+ 1) 

x(HTkWkyt  +  HTk+lWk+lyt+l)  (13) 

Substituting  (13)  into  (12)  and  simplifying  we  obtain: 

•*4+1  —  Xk  +  (HjWkHk  +  Hj+lW  k+i  H  k+i) 

x  Hj+lW k+i(yl+l  —  Hk+\xk)  (14) 

Eq.  (14)  is  of  the  form: 

Xk+i  =xk  +  Kk+i(y*k+l  -  Hk+lxk)  (15) 

where 

Kk+l  =  ( HjWkHk  +  HTk+lWk+lHk+lfl  HTk+lWk+ 1  (16) 


y(t)  =  h(x{t)) 


(19) 


where  W(t)  is  the  weight  on  the  process  noise  G(t )  and  u{t)  is 
the  input.  The  experimental  data  is  represented  by: 

yk=h(xk)  +  vk  (20) 

where  14  is  zero  mean  Gaussian  noise  associated  with  the  mea¬ 
surements,  whose  variance  is  a.  The  symbol  Xk  is  used  to 
approximate  the  actual  values  of  the  discretized  state  variables 
Xk .  This  approximation  is  essential  because  the  nonlinear  system 
does  not  have  an  exact  solution  like  the  linear  case.  The  objec¬ 
tive  then,  is  to  design  a  filter  quite  similar  to  Eq.  (14)  that  can 
predict  Xk+\  given  &k.  This  equation  is  of  the  form: 

34+1  =  i4  +  [K(xk)]M(yt  -  *(**))  (21) 


The  main  difference  between  Eqs.  (21)  and  (15)  is  that  the  gain 
matrix  (K(xk))  for  the  nonlinear  case  is  a  function  of  the  state 
variables  (xk)  and  hence  an  iterative  solution  procedure  is  essen¬ 
tial.  Upon  linearizing  the  model  equations  (Eqs.  (18)  and  (19)) 
an  equation  set  for  an  extended  Kalman  filter  can  be  obtained. 
The  linearized  form  of  the  matrices/  and  h  are  given  by: 


(Kfixity  u(t),  Q) 
dx(t) 


x(t)=xk 


(22) 


/mxm\ 

V  dx(0  J  x(t)=xk 


(23) 


For  known  values  of  G(t )  and  a  the  Riccati  equation  [41,45]  can 
be  used  for  the  filter  gain: 


K  _ _ Q  •  ^£+i^I+i _ 

k+i  ~  Q-Hk+lPk+lHj+l  +  Rk+l 
Pk+i  is  the  state  error  covariance  defined  by  [45]: 

Pk+ 1  =  £{++i+J+i ) 


(24) 


(25) 


Q  is  the  covariance  of  the  process  noise  (i.e.,  G(t)W(t)).  In  this 
work,  the  weight  matrix  W  is  held  constant  over  time,  for  sim¬ 
plicity  and  equals  the  identity  matrix  since  the  process  noise  is 
assumed  to  be  uncorrelated.  For  this  case,  the  covariance  of  the 
process  noise  is  defined  by  [45]: 

Q  =  £{WWt)  (26) 


Table  1 

List  of  parameters 
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Parameter 

Unit 

Anode 

Cathode 

Radius  of  the  particle  ( Rj ) 

m 

1 1  X  10'6 

12.5  x  10"6 

Rate  constant  (kj) 

moi-(1/2)  s-1 

4.907  x  10“5 

1.139  x  10“5 

Maximum  lithium  concentration  (LA1* ) 

mol  m~3 

51555 

30555 

Solid  phase  diffusion  coefficient  (D\j) 

2  -1 
m  s 

1.0  X  10“14 

3.9  x  10“14 

Mass  of  active  material  (mj) 

g 

15.92 

7.472 

Electrolyte  concentration  (co) 

mol  m~3 

1000 

Temperature  ( T) 

K 

298.15 

Faraday’s  constant  (F) 

Cmol-1 

96487.0 

Gas  constant  ( R ) 

J  mol-1  K-1 

8.314 

74+i  is  the  covariance  of  the  noise  in  the  measurement  (i.e.,  vk+i) 
given  by  [45]: 

Rk+i  =  E{vk+ivJ+l}  (27) 


The  Riccati  equation  is  a  result  of  a  transformation  that  seeks 
to  replace  a  differential  equation  with  an  equivalent  collection 
of  uncoupled  equations  of  lower  order  [43].  Two  limiting  cases 
are  considered  to  illustrate  the  capabilities  of  Eq.  (24):  case  (i) 
the  measurement  noise  is  much  greater  than  the  process  noise 
(i.e.,  74+i  Q )  and  case  (ii)  the  measurement  noise  is  much  less 
than  the  process  noise  (i.e.,  74+ 1  Q).  For  case  (i),  the  Riccati 
equation  (Eq.  (24))  simplifies  to: 


Kk+ 1  = 


Q  •  Pk+iHj+l 


Q  ■  Hu+\ Pk+\ rft+i  +  Rk+l 


0 


(28) 


In  other  words,  the  gain  matrix  is  not  updated  if  the  measurement 
noise  is  much  greater  than  the  process  noise.  For  case  (ii),  Eq. 
(24)  can  be  simplified  to: 


Q-Pt+ ig+i 

Q  ■  Hk+\ Pk+\ uj+l  +  Rk+ 1 


~  (Hk+\)~X 


(29) 


For  this  case,  the  gain  matrix  (Kk+ 1)  directly  maps  the  noise  (vk) 
in  the  output  (yk)  to  the  update  in  the  state  variable  (xk).  Thus,  the 
Riccati  equation  is  capable  of  filtering  the  measurement  noise 


out  while  retaining  the  process  noise.  The  choice  of  Q  and  74+ 1 
are  arbitrary — however,  determination  of  an  appropriate  value 
for  Q  from  the  variance  of  experimental  data  is  possible  [42,46] . 
More  details  on  these  tuning  parameters  are  provided  below. 


4.  The  single  particle  model 


The  single  particle  (SP)  approach  originally  developed  by 
Haran  et  al.  [31]  for  the  metal  hydride  battery  and  later  extended 
to  the  lithium  system  [32,33]  is  used  in  this  work.  This  model 
includes  the  assumption  that  each  electrode  of  a  lithium  ion  cell 
can  be  approximated  by  a  single  spherical  particle  whose  surface 
area  is  scaled  to  that  of  the  porous  electrode.  Further,  concentra¬ 
tion  and  potential  changes  in  the  solution  phase  are  ignored.  All 
parameters  are  held  constant  (see  Table  1)  and  thermal  effects 
are  assumed  to  be  negligible.  No  attempt  is  made  in  this  work 
to  model  capacity  fade  [30]  with  cycling  of  the  cell.  Most  of 
these  assumptions  can  be  readily  relaxed  [32]  at  the  cost  of  a 
more  sophisticated  model.  Fig.  1  provides  a  schematic  of  the  SP 
model  and  the  model  equations  are  presented  below. 

Diffusion  in  the  solid  phase  is  governed  by  the  Fick’s  laws 
[35]  written  in  spherical  coordinates: 


fcij  =  D 

dt  1  r~  dr 


(30) 


Fig.  1.  A  schematic  representation  of  the  single  particle  model. 
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chj 


t= 0 


=  c 


0 

1 J 


(31) 


( 


9£U 

dr 


=  0 

r= 0 


9cij\ 

9r  )  r=Rj 


ll 

nFDijaj 


(32) 


where  Jj  is  the  flux  entering  electrode  Eqs.  (1)— (3)  can  be 
volume  averaged  [36,37]  and  the  concentration  inside  the  solid 
phase  (c  \j)  can  be  expressed  in  terms  of  the  concentration  at  the 
surface  of  the  sphere  (ci  •)  and  the  average  concentration  inside 

the  sphere  (c?vf ).  The  resultant  equations  are: 

1  ’  J 


dcavg 

achj 

d  t 


+ 


15^U 


R 


(cTj  -  c\j)  =  0 


(33) 


/=o 


(34) 


Butler- Volmer  kinetics  [38]  is  employed  to  represent  the  charge 
transfer  reaction  at  the  surface  of  the  sphere: 


Fajkj 


max 


—  C 


exp 


—  exp 


0.5  F 

-gj-Wij  - 


(35) 


Finally,  since  none  of  the  above  variables  solved  for  (<p\j,  c\  ., 
c?vf )  are  directly  measurable  from  a  functioning  cell,  the  solid 

4  9  J 

phase  potentials  of  the  individual  electrodes  are  related  to  the 
cell  voltage  as  follows: 


T^cell  —  01,/?  01,  n  (36) 

The  open  circuit  potentials  (Uj)  are  empirical  functions  of  the 
solid  phase  concentration  (ci  ■  for  the  present  case)  obtained 
by  fitting  experimental  data  to  polynomial  expressions  [40]. 
The  expressions  presented  by  Ramadass  et  al.  [30]  are  used 
here: 


Uen  =  0.7222  +  0.13870,,  +  0.0290, °'5 


0.0172  0.0019 

-  + 


0 


n 


01.5 

A? 


+0.2808e(°'90-15e")  -  0.7984e(0'4465e"_0'4108) 


(37) 


5.  Online  estimation  of  the  state  of  charge 


The  SP  model  presented  in  the  previous  section  needs  further 
simplification  to  be  rewritten  in  the  form  of  Eqs.  (18)  and  (19). 
Eqs.  (33)  and  (34)  can  be  combined  and  the  resulting  equation 
can  be  rewritten  for  each  electrode  in  state  space  form  as  follows: 


r  dravg  1 
aclp 

- 1 

d  t 

dcavg 

UC1  ,n 

-  d  t  - 

1— 

3 

RpFctp 


3 

Rn  Fan 


(40) 


For  this  case  the  equivalent  of  Eq.  (19)  is  the  expression  for  the 
cell  potential  (VCell)  obtained  by  combining  Eqs.  (35)-(39)  to 
give: 


Veen  =  4>l,p  -  fa ,n  =  h(WSp(t),  Wlit),  Jp,  Jn)  (41) 

Vcell  corresponds  to  the  output  variable  y(t)  since  it  is  directly 
measurable.  The  vector  F  is  linear  for  this  model.  As  seen  from 
Eq.  (40),  we  have: 


3 

RpFctp 


3 

Rn  Fan 


(42) 


The  function  h  is  nonlinear,  and  is  dependent  on  the  state  vari¬ 
ables  (c?vg(t)  and  c?vg(t))  and  the  input  current.  The  experimental 
data  for  yceii  (represented  by  y%  in  Eq.  (20))  are  simulated  in  this 
paper  by  adding  a  zero  mean  Gaussian  noise  (vg)  with  a  variance 
of  0.1. 

The  state  of  charge  of  the  cell  (SOC)  is  represented  by  the  ratio 
between  the  average  concentration  of  lithium  (c?vf(t))  inside 
the  limiting  electrode  and  the  maximum  concentration  (c^x(t)) 
inside  that  electrode.  Most  commercial  lithium  ion  cells  are  lim¬ 
ited  by  the  performance  of  the  negative  electrode  [44] .  Hence, 
for  illustrative  purposes,  the  negative  electrode  is  assumed  to 
be  the  limiting  electrode  and  the  state  of  charge  based  on  the 
negative  electrode  ( 0n(t ))  is  given  by: 


(43) 


However,  the  technique  is  not  limited  to  negative  limited  elec¬ 
trodes;  an  equivalent  formulation  can  be  done  based  on  the 
average  concentration  of  lithium  inside  the  positive  electrode, 
for  the  case  of  positive  electrode  limited  cells. 


-4.656  +  88.6690^  -  401.1 1904  +  342.9090^  -  462.47108  +  433.4340^° 
-1.0  +  18.93302  -  79.53204  +  37.31106  _  73.0830^  +  95.960^° 


where 


(38) 


c 


1 J 


max 


C 


1 J 


(39) 


is  the  utility  coefficient.  At  the  end  of  complete  discharge,  0j 
reaches  a  value  of  1 .0  for  the  positive  electrode  and  0.0  for  the 
negative  electrode  in  an  ideal  cell. 


6.  Results  and  discussion 

The  process  noise  ( G(t)W)  is  a  measure  of  the  inaccuracy  of 
the  model,  while  the  measurement  noise  (vg)  is  a  measure  of 
external  disturbance  that  corrupts  the  actual  data  collected  from 
the  system.  Hence,  the  variance  of  the  process  noise  ( Q )  and  that 
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Fig.  2.  Cell  voltage  as  a  function  of  time  when  process  noise  is  negligible. 


of  the  measurement  (74+ 1)  can  be  used  as  tuning  parameters  of 
the  filter.  To  illustrate  this,  the  model  is  used  to  obtain  synthetic 
data  for  different  values  of  Q  and  74+ 1  •  The  synthetic  data  was 
generated  by  adding  zero  mean  Gaussian  noise  (G(t)),  with  a 
standard  deviation  of  10%,  to  the  discharge  curve  obtained  from 
the  single  particle  model.  Generating  data  by  adding  noise  to 
the  output  of  the  model  is  the  equivalent  of  having  an  error 
in  the  cell  voltage  due  to  faulty  measurements.  Hence,  a  good 
filter  should  eliminate  the  noise  due  to  the  measurements  and 
the  predicted  profile  for  the  cell  voltage  should  not  oscillate  as 
much  as  the  noisy  data.  Two  case  studies  are  considered  in  order 
to  illustrate  the  capabilities  of  the  present  approach.  In  the  first 
case,  all  the  noise  in  the  data  is  attributed  to  the  measurement 
error,  consistent  with  the  methodology  adopted  to  generate  the 
synthetic  data.  In  the  second  case,  an  incorrect  assumption  is 
made  that  all  the  noise  in  the  data  is  from  the  cell,  owing  to 
physical  changes  happening  within  the  system.  Fig.  2  shows 


the  estimated  cell  voltage  as  a  function  of  time  along  with  the 
synthetic  data  for  Case  1.  The  smooth  line  through  the  center 
is  the  predicted  profile  for  the  cell  voltage  from  the  cell.  This 
figure  illustrates  the  capability  of  the  EKF  approach  to  predict 
the  true  cell  voltage  filtering  the  noise  due  to  the  measurements 
out.  There  is  a  good  fit  between  the  data  and  the  predictions.  The 
tuning  parameters  used  to  generate  this  data  are  Q-  le— 6  and 
74+1  =  le— 2.  This  represents  the  case  where  the  process  noise 
is  insignificant  compared  to  the  measurement  noise.  This  is  true 
for  the  present  case,  since  all  the  disturbance  was  added  to  the 
true  potential  (i.e.,  the  actual  model  predictions)  as  measurement 
noise.  As  a  result  the  estimated  voltage  does  not  get  perturbed 
with  the  experimental  noise  seen  in  the  data.  This  verifies  the 
results  implied  by  Eq.  (28).  Also  shown  is  the  error  between  the 
two  curves  as  compared  to  the  bounds  set  by  the  experimental 
variance.  The  error  lies  well  within  the  bounds  of  the  3  —  a 
outliers  (shown  on  either  side  of  the  error  on  the  same  plot). 


O 

> 


Prediction  of  the  Cell  Voltage 


3.9 

3.8 

3.7 

3.6 


-  ,  1 — > - r 


I.  I 


s 


Q  =  1  e-3 
R  =1e-2 

k+1 


'V 


0 


500  1000  1500  2000  2500  3000  3500  4000  4500  5000 


Time  (sec) 


Fig.  3.  Cell  voltage  as  a  function  of  time  when  process  noise  is  considerable. 
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Fig.  4.  State  of  charge  as  a  function  of  time  when  process  noise  is  considerable. 


Fig.  3  shows  the  change  in  Vcen  with  time  for  Case  2.  For 
this  case,  we  assumed  that  all  the  noise  is  attributed  to  the  pro¬ 
cess  response,  i.e.,  the  inadequacy  of  the  model  to  reflect  all  the 
changes  that  happen  within  the  cell  and  hence  set  the  penalty  on 
the  model  ( Q )  to  be  le— 3,  which  is  considerably  higher  than  the 
previous  case.  The  value  of  74+ 1  is  the  same  as  in  the  previous 
case.  These  values  place  less  emphasis  on  the  experimental  noise 
and  attribute  the  oscillations  in  the  data  to  unmodeled  dynamics. 
For  this  case,  one  would  require  that  the  predicted  cell  voltage 
follows  the  changes  in  experimental  Vceii  more  accurately,  since 
the  oscillations  in  this  case  are  assumed  not  to  originate  from 
measurement  error.  As  seen  from  Fig.  3  the  predicted  profile 
for  the  cell  voltage  closely  follows  the  changes  in  the  data  as 
predicted  by  Eq.  (29). 

For  estimation  of  the  state  of  charge,  Eqs.  (33)  and  (34)  were 
used  along  with  the  other  model  equations  to  obtain  the  SOC 


inside  the  limiting  electrode  (0n(t))  as  shown  in  Eq.  (43).  Fig.  4 
shows  the  state  of  charge  calculated  using  Eq.  (43)  as  a  function 
of  time  for  the  second  case  (i.e.,  Q=  le— 3  and  74+i  =0.01). 
This  response  is  noisy  and  the  error  is  about  ±10%  as  shown 
by  the  outliers.  The  high  error  in  the  state  of  charge  estimates 
reflects  our  incorrect  assumption  for  Case  2  that  the  oscillations 
in  the  data  are  due  to  physical  changes  into  the  cell.  Note  that 
the  synthetic  data  were  generated  by  adding  measurement  noise 
to  the  true  cell  potential.  Hence,  this  result  is  expected.  In  other 
words,  no  process  noise  was  used  in  generating  the  synthetic 
data  and  all  oscillations  in  the  Vceii  versus  t  data  are  because 
of  the  measurement  noise  added  to  the  model.  Since  the  tuning 
parameters  used  in  Case  2  penalize  the  model  heavily  despite  its 
reasonable  accuracy,  the  inaccuracy  due  to  an  error  in  judgment 
is  reflected  in  the  SOC  predictions.  This  inference  is  further 
verified  in  Fig.  5,  where  the  process  noise  is  set  to  a  negligible 


Estimation  of  the  State  of  Charge  Using  Q  =  1e-6 


Fig.  5.  State  of  charge  as  a  function  of  time  when  process  noise  is  negligible. 
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Fig.  6.  Comparison  of  estimates  to  experimental  data. 


value  (Case  1)  and  the  scatter  in  the  response  is  attributed  to  the 
measurement  noise.  The  state  of  charge  estimate  for  this  case  is 
much  smoother  and  follows  closely  the  values  used  in  the  model. 
The  error  in  the  state  of  charge  predictions  is  about  ±2%. 

In  order  to  verify  the  validity  of  the  model,  experimental 
data  from  Sony  18650  Cells  are  compared  to  the  estimated  cell 
voltage  and  state  of  charge  in  Fig.  6.  The  discharge  curve  was 
obtained  from  at  25  °C  in  an  Arbin  battery  testing  system  (Col¬ 
lege  Station,  TX)  at  a  discharge  current  of  1  A,  after  the  cell 
was  charged  to  its  full  capacity  of  1.78  Ah.  The  data  obtained 
from  the  experimental  discharge  curve  were  used  to  update  the 
model  estimates  at  each  instant  of  time.  The  experimental  state 
of  charge  at  each  instant  of  time  was  obtained  by  integrating 
the  charge  removed  from  the  cell.  The  values  for  Q  and  74+ 1 
were  held  at  le— 6  and  0.01,  respectively.  Fig.  6  shows  very 
good  agreement  between  the  experimental  cell  voltage  and  the 
predictions. 

7.  Advantages  of  the  current  approach 

The  methodology  outlined  in  this  work  does  not  involve 
recursive  nonlinear  regression.  The  conventional  estimation  pro¬ 
tocols  invoke  a  nonlinear  regression  technique  [39,42],  which  is 
a  complex  iterative  process.  In  order  to  obtain  V  parameter  val¬ 
ues,  one  needs  to  solve  a  minimum  of  V  equations — and  since 
the  Kalman  filter  approach  for  the  estimation  of  V  parameters 
involves  inverting  a  matrix  of  size  no  bigger  than  ‘n  ,  this  is  the 
most  optimal  algorithm  to  estimate  parameters  online.  Whereas 
the  competitive  non-linear  least  squares  algorithm  requires  about 
2.8  s  of  computational  time  on  a  Pentium  IV,  2  GHz  processor,  in 


order  to  obtain  the  state  of  charge  at  a  particular  time,  the  Kalman 
filter  approach  presented  here  takes  about  4  ms  for  an  equiva¬ 
lent  estimate.  This  rapid  estimate  of  the  state  of  charge  enables 
one  to  follow  closely  the  performance  of  the  cell  as  shown  in 
Figs.  2-5.  Another  advantage  of  the  current  approach  is  that  it  is 
simple  enough  to  be  implemented  in  a  real  time  device,  since  it 
spares  the  tedium  of  solving  a  complex  nonlinear  optimization 
problem. 

Incorporation  of  a  physics  based  model  along  with  the 
Kalman  filter  based  estimation  technique  improves  the  predic¬ 
tive  capability  of  the  model,  in  addition  to  providing  a  physical 
insight  to  the  state  of  charge  estimates.  For  example,  the  state  of 
charge  in  the  present  approach  is  the  ratio  between  the  average 
concentration  to  the  maximum  concentration,  of  the  active  mate¬ 
rial  within  the  electrode  that  limits  the  performance  of  the  cell, 
unlike  the  empirical  values  predicted  by  the  models  in  the  liter¬ 
ature.  Plett  [19-21]  introduced  some  models  to  predict  the  state 
of  charge  of  a  lithium  ion  cell  accurately  using  the  extended 
Kalman  filter  approach.  However,  results  from  his  work  [20] 
show  that  it  is  essential  to  add  an  arbitrary  number  of  poles 
and  zeros  to  the  transfer  function  used  to  represent  the  cell  in 
the  model,  in  order  to  obtain  a  very  good  fit.  This  approach  is 
as  empirical  as  the  circuit  based  approach  proposed  by  others 
[14-16].  The  state  of  charge  predicted  by  the  empirical  mod¬ 
els  does  not  have  a  direct  correlation  with  the  dynamics  of  the 
cell  as  do  the  predictions  from  rigorous  models  based  on  the 
transport  and  kinetic  equations.  In  addition  to  this,  results  from 
Fig.  6  show  that  a  physics  based  models  can  be  used  to  obtain 
estimates  of  the  state  of  charge.  A  Monte  Carlo  simulation  was 
performed — generating  100  sets  of  synthetic  data  and  the  max- 
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Trial 

Fig.  7.  Maximum  error  between  SOC  from  model  and  synthetic  data  for  100 
trials. 

imum  error  in  the  estimation  of  the  state  of  charge  for  each  run 
is  plotted  in  Fig.  7.  The  results  show  that  for  all  the  cases,  the 
estimates  lie  well  within  ±2%. 

8.  Conclusion 

A  model-based  extended  Kalman  filter  was  developed  based 
on  a  simple  electrochemical  model  of  a  lithium  ion  cell.  The 
state  of  charge  of  the  cell  is  monitored  continuously  as  a  func¬ 
tion  of  time.  The  predictions  were  compared  with  synthetic  data. 
The  Kalman  filter  approach  represents  the  optimal  algorithm  for 
recursive  identification  [16].  The  SP  model  adds  the  information 
about  the  time  evolution  of  the  parameters  characterizing  the 
system  state.  Hence,  the  approach  presented  here  provides  esti¬ 
mates  of  the  SOC  of  the  cell,  that  are  as  good  as  any  empirical 
methodology  adopted  in  the  literature,  while  providing  a  good 
physical  insight  into  the  system.  To  summarize,  the  methodology 
outlined  here  greatly  reduces  the  computational  time  involved  in 
obtaining  the  state  of  charge  from  a  first  principles  based  model, 
while  adding  physical  insight  into  otherwise  equivalent  circuit 
based  models  implementing  the  filtering  technique. 
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